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Редукція задачі Алексідзе до рівняння сили тяжіння 
(Представлено академіком НАН України В.І. Старостенком) 

З метою аналітичного продовження в глобальних областях значень сили тяжіння, які є значеннями модуля 
градієнта потенціалу сили тяжіння (МГІПСТ), сформульована у вигляді інтегрального рівняння нелінійна 
гранична задача Алексідзе для рівняння Лапласа. Для її розв'язання у вигляді потенціалу простого шару на 


поверхні Ляпунова слід визначити невідому густину простого шару. На основі аналітичних властивостей функції 
МГПСТ задачу Алексідзе зведено до нелінійного інтегрального рівняння сили тяжіння. 


Методи трансформації потенціальних полів дають змогу розв'язувати з достатньою точністю 
будь-яку задачу за умови задання вхідної інформації в локальних областях певної малої міри. На- 
явні дані гравітаційних і магнітних спостережень як різниці приростів модуля градієнта потенціал- 
у сили тяжіння (МГПСТ) можна використати при розробці схем трансформацій потенціалу в гло- 
бальній області. 

Для глобальних побудов густинних моделей земної кори, які виконують на основі даних регі- 
ональних спостережень, потрібно розв'язати задачу аналітичного продовження аномалій сили тя- 
жіння (1 |. В праці |2) виокремлено два альтернативні напрямки, що базуються на аналізі характе- 
ристичних властивостей МГПСТ. Один з них зводиться до розв'язання зовнішньої задачі Діріхле 
для лінійного диференціального рівняння типу Гельмгольца, яке моделює значення МГПСТ в об- 
ласті, не зайнятій тяжіючими масами, з невідомим змінним коефіцієнтом (|З). 

Перехід до задачі відновлення потенціалу усунув необхідність обчислення наступних набли- 
жень як коефіцієнтів рівняння сили тяжіння, так і самих значень сили тяжіння, бо останні тепер 
можна знайти не лише з розв'язання задачі Діріхле для рівняння сили тяжіння, а і з безпосередньо- 
го диференціювання відновленого потенціалу. Через це задача відновлення потенціалу за значен- 
нями МГПСТ, яка не належить до класичних задач гравіметрії, набуває особливого значення. 
Один з можливих способів її вирішення розроблений в |4|, а інший - гранична задача Алексідзе 
для рівняння Лапласа - викладено в (5, 61. 


Постановка задачі. Нехай у - обмежена область точок тривимірного евклідового простору, 
зайнята масами Землі, у" - необмежене доповнення до у , що вільне від тяжіючих об'єктів, ду - 
границя множин у і у" , що ототожнена з фізичною поверхнею Землі. В прямокутній декартовій 
системі координат Ох;,хь,х, з початком у центрі Землі, осі Ох, , Ох, лежать в екваторіальній площи- 
ні, авісь Ох, співпадає з віссю її обертання. Точки області у позначмо через 2 - С ом 1 У, 
доповнення у - через х- (ззоьве у".Маси М (2), є є у" згустиною АМ (2) Е са всереди- 
ні Землі генерують у просторі у" поле сили тяжіння з потенціалом 

орел р де | хеу 
а Ра 4 0,  хеу" 
де 7 - гравітаційна стала, Ох) -0.5а (х Я хо) - потенціал відцентрової сили, шо - модуль вектора 


кутової швидкості Землі. Напруженість гравітаційного поля (значення ПГПСТ згідно з |2|) дорів- 
нює: 


віч рутіо) « У ОО дк о я) ау 


дх, (х) : ее і а : 
де Р б соб(п, Ху , К - 1,2,3 - напрямні косинуси одиничного вектора п(х) внутрішньої нормалі 
п 
до еквіпотенціальної поверхні «И/(х): М/(у) - Сх, яка проходить через точку х. 
Для такої моделі середовища необхідно вирішити нелінійну граничну задачу Алексідзе для 


рівняння Лапласа: знайти функцію (2), х є у", яка задовольняє всередині необмеженої замкну- 


тої області У з у" ду рівнянню Лапласа ЛИ/(х)- 0, х є у", якщо в будь-якій точці ляпуновсь- 
кої границі ду області і в нескінченно віддаленій точці вона задовольняє умовам: 
2 
зро) а а? (х), х єдбу, М (х)-з 0 при Їх -»со, (2) 
ст бХ, 
де «(х) - задана неперервна функція вигляду (1). 
Гармонічну в області у" функцію й | хє "й дня як потенціал простого шару 


б 


я 5, 


з невідомою густиною с(х), хеду РА, хоча загалом може бути більш гладкою функці- 


5р хе" (3) 





єю), поширеного на поверхні ду Ляпунова. Щоб вивести рівняння, з якого відновлюють невідому 
густину с(х) за заданими на поверхні бу значеннями «(х) МГПСТ та для обчислення значень 
«(х) в точках області у" при знайденій густині, визначмо, виходячи із зображення (3) потенціалу, 
квадрат модуля його іп 

1 
)- вгтади (х | з 
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ра8 сп) я Хровк Хо Пе Пк 45 а8.. 
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Цей вираз стане компактнішим, якщо ввести в обіг одиничні вектори р і 4, спрямовані з точ- 








ки х в точки С і 7 відповідно, що пробігають по поверхні ду. Їх складові матимуть вигляд 




















ду ду 


р, з сов(рух,)- ра , 4; 2 сов(х,,4)- ет і 
о езвводнь оон КЕ -Уоюні і не «У ба і б б | 
З урахуванням цих позначень вираз для є 2) б о 
золот з імаеь з со5(р, 45, аб, (4) 





Проаналізуємо аналітичні особливості функції б модуля градієнта потенціалу простого ша- 
ру (4). Першою з них є неперервність, що випливає з наступної леми. 

Лема 1. Модуль градієнта потенціалу простого шару є неперервною функцією точки хєду, 
яка рухається по поверхні ду Ляпунова. 

Доведення. Слід переконатися в існуванні інтеграла (4), поширеного на частину 
дух) збупг5 (2 є) поверхні ду всередині кулі 5 ее) Ляпунова радіуса є » 0 з центром в точці 
х є ду, оскільки його існування в іншій частині Оу ду(х) не підлягає сумніву. Доведення існував- 
ня інтегралу по ду(х) полягає у відшуканні нерівностей, що витікають з локальних властивостей 
поверхонь Ляпунова, які властиві й поверхні ОУ. 

Лема 2. Квадрат модуля градієнта потенціалу простого шару, поширеного на сфері радіусом 
р з одиничною поверхневою густиною с(х) -1, хе ду, становить 








»р 
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Доведення. Розглянемо спочатку випадок, коли точка розташована зовні сфери бу радіусом 





сов р.) с о У МВА. 
ту» Рене ЗМР В о, 








Р . Покладемо в точці х початок допоміжної системи координат ОС,2,65, спрямувавши вісь ОС, 
крізь центр сфери. Тоді центр сфери матиме координати (0,0, а), де а -|х - апліката центра 


сфери, а її рівняння - б Є; в (2, б а) - р". В полярних координатах рівняння набуває вигляду: 


г (с, 4)- 2аг(ф, А)совфча? - р'-б, 
оскільки Є з Мо, 4)віпфсовл, фе го, А)віпфзіпа, бе по, Л)созф. Звідси випливає, що о, 2) -лд), 
причому "(ф) засозфІі ра -а? віп Ф,0«ФеФ,, де Ф, з агсвіп р/а. Виходячи з цього виразу, 
рівняння сфери можна однозначно поділити відносно напряму осі ОС, на дві частини - ,,верхню" 1 


»нижню?, що задаються, відповідно, рівняннями: 


г. (р) - асовф чн р? - авіп' ф, 0«ФАФ,, 
кр) - асовф-/р' - а? віп'ф, 0сфФаФ,. 


Знайдемо найменший елемент 45 сфери бу, параметризованої полярними координатами, для 





де, да, дб, 
чого скла др др ді заїї й щі 
демо матрицю і за її елементами підрахуємо коефіцієнти Гауса сфери 
де, 05, дб; 
0 ОА дА 


ду , отримавши при цьому 


45 - к(ф, А) Ї"(ф)-- Р зіп дааїЛ.. 


Звідси з урахуванням зображень ,,верхньої" і ,нижньої" частин сфери, отримуємо значення 
елементів поверхні ,,верхньої" 1 ,,нижньої" частини сфери відповідно, у вигляді 











2 : 2 5 
48, а (ф)рвіпф з крабів (ф)рзіпф ЯРУ 
р -а"звіпф р'-а"звіпф 
Враховуючи, що Я - "ф) та виразивши спрямовуючі косинуси одиничних векторів р і 4 че- 


2 
Фо 





. 4 
р дане ре я , що й треба було довести. 
она 


Нехай точка х лежить на сфері бу. Покладемо в ній початок допоміжної сферичної системи 


рез ф і 4, отримаємо я (х)- 


координат та отримаємо рівняння сфери "(ф) е2рсо8ф,0«фФХ«л/2,1, відповідно, елемент повер- 


хні д5 « МЕС - К'адад « Ар' віпфсозфафай,, де Е -« Е(о, 4), С - С(ф, А), Е « Е(ф,Л) - коефіціє- 


нти Гауса поверхні ду. Підставляючи останній вираз у вираз для 8" (х) , одразу отримуємо 
2 


2 1 Р 1 
(чі Гдеевняння 2, 
417 г" (ф) 4 


що доводить лему для випадку, коли х є бу. 
Якщо точка х розташована всередині сфери бу, покладемо початок допоміжної сферичної 


системи координат в точку х з полярною віссю, що проходить крізь центр сфери, та отримаємо 
рівняння сфери для точок її , верхньої" і нижньої" частин 


г (ф)  асовф-нр' -а"'звіп'ф, 0 «фХл/2, 
г (д) асовф-Ар' -а"'віп'ф, -лі2«ф«0, 
деасз Їж - апліката центра сфери. Звідси після нескладних перетворень цих зображень, отримаємо 
а"(х)- 0 при | «фіз 
Функція а "(х) (4) є розривною і має розрив неперервності при перетині точкою х поверхні 
ду . Щоб знайти величину розриву, позначмо через 85 (х) пряме значення квадрата МГПСТ, якщо 


точка х лежить на поверхні ду, а через «Цх.) і а.) - граничні значення ач); якщо точка х 


прямує до поверхні ду поза та зсередини області у по нормалі до ду. Величину розриву роз- 
криває теорема. 

Теорема 1. Нехай простий шар поширений на поверхні Ляпунова бу, що охоплює область у . 
Якщо густина ох), х єду потенціалу простого шару неперервна, то квадрат МГПСТ а "(х.) при 
наближенні точки х, в напрямі нормалі п(х) до точки х на поверхні ду ззовні та зсередини об- 
ласті у має різні границі: 

іти" (ку) а р а, го) 
4 4ло 3 о 4| 


Хо-Р?хеду ду 





Хоєпбу 
Р обро а с(х) соз(п, р) 45 2 
шиє (хо) 28; (х)- С (2) Ї 2 с(Є) РОН ЕВ 
4 4л ду ре -4| 
Хо ?хеду 
хоєпоу 
Доведення. Щоб довести теорему, застосуємо означення МГПСТ та твердження з публікації 
131: якщо густина потенціалу простого шару (3) є диференційовною, граничні значення частинних 
похідних потенціалу 1-го та 2-го порядку визначаються величинами 
Ох.) 1 сх, -б, 


і 1 
дх, Ь 4я3 ре б Є сім, 5 5 х)совіх, й п) у 











си 1 с3бх, -6,) -ре-еЇ іогобоі 
ШОВ о НМ обом, робо 


ОТ) іму 2 З -єДх, -6) 


дх, дх, я а. 
і бу 








с(2) с(х))д8. ВРЮ із), і,/ «12,3, 
Хі 
які можна отримати з подання зозрни И(х) у вигляді 


нед Очеамя, ков 


бу 


а8. 
хоб 


якщо замінити безпосереднє диференціювання за х, диференціюванням по внутрішній нормалі п 





2 








до поверхні ду в точці х та врахувати тотожність 


і, хеу" 
| рел суч, 2112, хеду, 
4л з дп с -4 4 бо 


що випливає з тотожності Гріна при У(х) ее 


Дійсно, при хє п/п у" маємо 


йт (асо) кб) ев б - 


хо-»хеду га дху ДЛ бу дп Р і 4 





«бобу ре оду, і), 





і аналогічно при хє пу маємо 


інь У Й о ан Па) 0 ре ау 


Хо?хебу ОХ 4 Дя - РАЗУ 





Ка! 
1 г с(2) ос(тп) 
ноя зебії 
де 86(2)- бл" асо ог зон 











Звідси при с(х)е 1, ва ню. 


хєду та з умови для сфери - 45, --- отримуємо 
47 5, х-4| 2 


2(2)- єс (х)- 3/4, х є ду, що узгоджується з лемою 2. 


Оскільки функція а" (х) при наближенні точки х до поверхні ду має розрив (доведений тео- 
ремою 1), а при русі тієї ж точки по поверхні бу вона неперервна (як стверджує лема 1), то для ви- 
конання граничних умов (2) задачі Алексідзе необхідно, щоб виконувалась рівність 
2: (х) - а" (х), х є бу , тобто задані на границі ду значення сили тяжіння дорівнюють значенням в 
зовнішньому просторі. 

Зауважуючи зображення для є: (х) в теоремі 1, отримуємо щодо невідомої густини потенціалу 
простого шару (3) нелінійне інтегральне рівняння сили тяжіння: 

аа (х)ч сх) Ї совіль д) с(2)45. - | 5 Г оз) с). со5і р, 4)45, 45; зах) хеду. (5) 

4 4л 5, жо 16л х-й ж - п 


Розв'язок с(х), хєду інтегрального рівняння (5) еквівалентний розв'язку задачі Алексідзе з гра- 








дуду 





ничними даними на поверхні ду Ляпунова, оскільки за будь-якого вибору густини потенціал (3) 


простого шару задовольняє в області у" рівнянню Лапласа, а знайдене з рівняння (5) значення 
густини забезпечує виконання граничної умови (2). Питання розв'язності задачі Алексідзе по суті 
редукується до вияснення умов існування, єдності та стійкості розв'язку нелінійного рівняння (35). 
Це рівняння можна спростити, виходячи з того, що (наслідок леми 2) внутрішня точка не тяжіє 
до однорідного сферичного шару. Отже, можна припустити, що існують не лише для сферичних, 
але й для інших шарів такі розподіли густини, які не притягують внутрішніх точок. У справед- 
ливості припущення переконаємось, якщо проаналізуємо взаємодію внутрішніх точок області з 
простим шаром в точках границі бу області, що є поверхнею Ляпунова, з густиною, яка задоволь- 
няє інтегральному рівнянню 
іо р ру, нах) 0, хеду, (6) 
4 4ло 5, | хо 4 
яке виходить з (5) з урахуванням виразу (4). 
Якщо густина с(х) простого шару є розв'язком рівняння (6), то цей шар створює потенціал 


(ох) , Значення якого не змінюються всередині області у . Дійсно, умова (6) відповідає граничній 
. й й б є га 
рівності 2.(х)- 0, х є ду, яка, в свою чергу, еквівалентна співвідношенню (тайм (2) -0, хєду. 


Звідси випливає, що гармонічна функція И(х), хє у" в кожній точці границі ду Ляпунова задо- 

дих) | М Че - 

вольняє умову ------ 0, Кк 12,3, тобто потенціал Ух) з соп5і сталий у всій області хєу . 
Ху 


Через це внутрішні точки області у не тяжіють до простого шару з густиною, яка задовольняє 





умові (6). 
Тепер додаймо вирази (5) 1 (6), в результаті чого отримаємо нове інтегральне рівняння 
1 гр) сп) 
с"(х)ч со5| р, 445,45. «227 х) хєду, (7) 
( ) Дл? Па РАВ ( М. 7 б ( ) 





розв'язок якого зумовлює розподіл сили тяжіння «(х) в необмеженій області у , що відповідає 
постійному значенню потенціалу (ох) всередині області у . Розв'язки рівнянь (5) 1 (7) допомага- 
ють визначати не лише потенціал (2), х є у", а й значення МГПСТ в будь-якій точці необме- 


женої області у". Їх можна обчислити як з (4), так і за значно зручнішою формулою 
2 


Рани о у Бл ЗРРУО 


3 
2 3 
16 ка ду х- 4 





Зведення задачі Алексідзе з граничними даними (2) на поверхні Ляпунова до розв'язання не- 
лінійного інтегрального рівняння сили тяжіння (5) дозволяє вивчити питання розв'язності та 
єдиності її розв'язків, ефективно знаходити числові наближення для областей складної форми. 
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Ведисіїоп ої пе АЇехідге ргобіет їо Пе сгауїсу едпабіоп 


Рог ригрозез оГ ап апабуїіс сопіїпиатоп о єтауйу уаїез аз а тоашіе оГ єтауйу роіепій! єтадіепіз (МСРСОІ іп яіоРаї! 


агеаз а попіїпеак Боипаагу-уаїе Аехідге ргобіет /ог Ше ІГаріасе'5 едиапоп із /огтиатей. То 50Їуе ії аз а 5ітріє 
Ідуек роіепііаі оп Те Гуарипоу 5 хинасе й із песеззагу 1о ае/їпе Ше ипКпомтп 5ітріе Іауег аеп5іу уаїнез. Оп Ше Базі 
оГ апаЇутіса! /еаїшгез о МСРОС /ипстіоп, Пе Аехіаге ргобіет із гедиседй іо Ше 50Ішіоп ої а попіїпеаг іпіеята! ятауйу 
едиаїоп. 
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